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ABSTRACT 



^"^ ' Aims. We present aperture- synthesis imaging of the red supergiant Antares {a Sco) in the CO first overtone lines. Our goal is to probe 

^^ , the structure and dynamics of the outer atmosphere. 

Methods. Antares was observed between 2.28 jim and 2.31 jim with VLTI/AMBER with spectral resolutions of up to 12000 and 
angular resolutions as high as 7.2 mas at two epochs with a time interval of one year 

Results. The reconstructed images in individual CO lines reveal that the star appears diff'erently in the blue wing, line center, and red 

C/j \ wing. In 2009, the images in the line center and red wing show an asymmetrically extended component, while the image in the blue 

^ . wing shows little trace of it. In 2010, however, the extended component appears in the line center and blue wing, and the image in the 

f^ ■ red wing shows only a weak signature of the extended component. Our modeling of these AMBER data suggests that there is an outer 

I ' atmosphere (MOLsphere) extending to 1.2-1.4 R^, with CO column densities of (0.5-1) X 10^^ cm~^ and a temperature of ~2000 K. 

O , The CO line images observed in 2009 can be explained by a model in which a large patch or clump of CO gas is infalling at only 

^ . 0-5 km s~\ while the CO gas in the remaining region is moving outward much faster at 20-30 km s~^ The images observed in 2010 

(/^ ■ suggest that a large clump of CO gas is moving outward at 0-5 km s~\ while the CO gas in the remaining region is infalling much 

C^ ' faster at 20-30 km s~^ In contrast to the images in the CO lines, the AMBER data in the continuum show only a slight deviation 

from limb-darkened disks and only marginal time variations. We derive a limb-darkened disk diameter of 37.38 + 0.06 mas and a 

power-law-type limb-darkening parameter of (8.7 + 1.6) x 10"^ (2009) and 37.31 + 0.09 mas and (1.5 + 0.2) x 10"^ (2010). We also 

obtain an eff'ective temperature of 3660 + 120 K (the error includes the eff'ects of the temporal flux variation that is assumed to be the 

same as Betelgeuse) and a luminosity of log L^/L© = 4.88 + 0.23. Comparison with theoretical evolutionary tracks suggests a mass 

of 15 + 5 M© with an age of 1 1-15 Myr, which is consistent with the recently estimated age for the Upper Scorpius OB association. 

Conclusions. The properties of the outer atmosphere of Antares are similar to those of another well- studied red supergiant, 

Betelgeuse. The density of the extended outer atmosphere of Antares and Betelgeuse is higher than predicted by the current 3-D 

convection simulations by at least six orders of magnitude, implying that convection alone cannot explain the formation of the ex- 

'^- \ tended outer atmosphere. 

^*^ I Key words, infrared: stars - techniques: interferometric - stars: supergiants - stars: late-type - stars: atmospheres - stars: individual: 

^ * Antares 
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1 ■ Introduction For understanding the wind acceleration mechanism, it is im- 

^ portant to study the region between the photosphere and the in- 

(^ The mass loss in the red supergiant (RSG) stage significantly af- nermost circumstellar envelope, where the stellar winds are ex- 
■ ■ ' fects the final fate of massive stars. For example, the mass loss pected to be accelerated. Various observations suggest the com- 
in the RSG stage seems to be a key to understanding the progen- plicated nature of this region. Tuthill et al. a997 ) have discov- 
itors of the most comn ion co re-collapse supemovae (SNe) Type ered asymmetric structures on the surface of three well-studied 
Up (e.g., Smartt et al. 120091) . However, the mass loss in RSGs RSGs, Betelgeuse (a Ori), Antares (a Sco), and a Her. They 
is still poorly understood. As Harper (2010l stresses, there are interpret these asymmetric structures as the presence of one to 
currently no satisfactory theories for the RSG mass loss. Even three hot spots. The recent near- IR inte rferometric observations 
the driving force of the stellar winds has not yet been identified, of Betelgeuse by Haubois et al. (I2009D also show surface inho- 

mogeneities. The detection of H2O in Betelgeuse and Antares at 
12 jjm (Jennings & Sada 1998), as well as in the near-IR and in 



Send offprint requests to: K. Ohnaka the 6 yum region (Tsuji 2000ai 2000b), suggests that there is an 
* Based on AMBER observations made with the Very Large outer atmosphere extending to 1.3-2 stellar radii with tempera- 
Telescope Interferometer of the European Southern Observatory ^^j-^s of 1500-2000 K, the so-called MOLsphere as coined by 
Program ID: 083.D-0333(A/B) (AMBER Guaranteed Time j -- /l oQQQbl ) 
Observation), 085.D-0085(A/B) ^ ^ ^' 
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The radio observations of Betelgeuse and Antares also show 
neutral gas extending to several stellar radii (Lim et al. 19981 
Harper et al. 2010). On the other hand, spatially resolved ob- 
servations of RSGs in the UV, as well as in the Ha line, sug- 
gest that the chromosphere also extends to several stellar radii 
(Gilliland & Dupree [T9961 White et al. [T9821 Hebden et al. 
1987). Harper et al. (2001) propose that the chromosphere has 
a small filling factor embedded in the neutral/molecular gas. 
However, the physical mechanism responsible for this inho- 
mogeneous, multicomponent outer atmosphere is by no means 
clear. Inhomogeneous structures are also found farther out from 
the star. The near- and mid-IR imaging of Betelgeuse (K ervella 
et al. 20091 [2QT1 ) and Antares (Bloemhof et al. [T9g4l [T9951 
Marsh et al. 1200 II) reveals a clumpy structure of the circumstellar 
envelope. 

For a better understanding of the mass-loss mechanism in 
RSGs, it is essential to probe the dynamics of the inhomoge- 
neous, multicomponent outer atmosphere. We carried out high 
spatial and high spectral resolution observations of Betelgeuse 
in the CO first overtone lines near 2.3 yum using VLTI/ AMBER 
(Ohnaka et al. 2009 ; 201 1 , hereafter Paper I and II, respectively). 
With the stellar disk spatially resolved with an angular resolution 
of 9.8 mas (the highest resolution ever achieved at any wave- 
length for this well-studied star), AMBER's high spectral reso- 
lution of up to 12000 enabled us to spatially resolve, for the first 
time, the gas motions over the surface of a star other than the 
Sun. These AMBER observations have revealed that the star ap- 
pears diflTerently in the blue and red wings of the individual CO 
lines due to an inhomogeneous velocity field in the photosphere 
and MOLsphere. 

In this paper, we present high spectral and high spatial reso- 
lution VLTI/AMBER observations of another well-studied RSG, 
Antares (M1.5Iab-b), in the CO first overtone lines. Antares is 
similar to Betelgeuse (M2Iab:) in terms of the spectral type, 
although it is slightly less luminous and less massive than 
Betelgeuse (see Sect. 13.21) . The mass-loss rate of Antares is esti- 
mated to be ~2 X 10"^ MoPERYR (Braun et al. [2QT2b , which is 
com parable to that of Betelgeuse (~3 x 10"^ MoPERYR, Harper 
et al. 1200 lb . However, the IR excess and the 10 yum silicate emis- 
sion of Antares are much weaker than those of Betelgeuse (e.g., 
Monnier et al. 119981 Verhoelst et al. 120091) . Our goal is to spa- 
tially resolve the gas motions in the outer atmosphere in this 
second RSG and to examine whether the findings for Betelgeuse 
are common among RSGs. The paper is organized as follows. 
The observations and data reduction as well as image recon- 
struction from the interferometric data are described in Sect.[2l 
We present the results in Sect. [3] followed by the modeling of 
the data in Sect. |4] and the interpretation of the results in Sect. [3 
Conclusions are given in Sect. [6l 



2, Observations 

2.1. AMBER observations and data reduction 



AMBER (Petrov et al. 12007b is a VLT spectro-interferometric 
instrument that operates in the near-IR (1.3 — 2.4 yum) and com- 
bines three 8.2 m Unit Telescopes (UTs) or 1.8 m Auxiliary 
Telescopes (ATs). AMBER records spectrally-dispersed inter- 
ferograms with a spectral resolution of up to 12000. AMBER 
measures the amplitude of the Fourier transform (so-called vis- 
ibility or visibility amplitude), closure phase (CP) and diflTeren- 
tial phase (DP), with the usual (i.e., spatially unresolved) spec- 
trum. The CP is the sum of the measured Fourier phases around 
a closed triangle of baselines (i.e., (fu + (f23 + ^31), which is not 



aflTected by the atmospheric turbulence. Therefore, it is essential 
for aperture- synthesis imaging in optical/IR interferometry. The 
CP is zero or n for point- symmetric objects and non-zero and 
non-TT CPs indicate an asymmetry of the object. The DP provides 
information about the wavelength dependence of the photocen- 
ter shift of the object within spectral features. 

Antares was observed on 2009 April 24 and 25 and, one year 
later, on 2010 April 10 and 11, with AMBER using the EO-GO- 
HO (16-32-48 m) array configuration and the EO-GO-Il (16- 
57-69 m) configuration (Program IDs: 083.D-0333, 085.D-085, 
P.I.: K. Ohnaka). As in Papers I and II, we used the high res- 
olution mode in the /sT-band (HR_K) with a spectral resolution 
of 12000 and covered wavelengths from 2.28 to 2.31 yum to ob- 
serve the CO first overtone lines near the (2,0) band head. We 
detected fringes on all three baselines without using the VLTI 
fringe tracker FINITO (Antares saturates FINITO) thanks to the 
extremely high brightness of Antares (K = -4.1). In 2009 and 
2010, we observed on two half nights each and obtained 65 and 
44 data sets, respectively. The data sets taken more than ~2 min- 
utes apart were treated as separate data sets, because Antares 
is strongly resolved, and the visibility varies noticeably even 
with a slight change in the projected baseline length (see Fig. [3]). 
Figured] shows the uv coverage of our observations in 2009 and 
2010. All data sets were taken with a Detector Integration Time 
(DIT) of 130 ms. A summary of the observations is given in 
Tables|Al]and|A21 

For the reduction of the AMBER data, we used amdlib ver. 
2.2B which is based on the P2VM algorithm (Tatulli et al. ,2007J . 
In order to improve the signal-to-noise ratios (S/N), we applied 
a binning in the spectral direction to all the raw data (object, 
dark, sky, and P2VM calibration data) using a running box car 
function as described in Paper I. The binning down to a spectral 
resolution of 8000 provided an S/N suflftcient for the image re- 
construction. One of the parameters in the reduction with amdlib 
is the frame selection criterion. In each data set, we checked for a 
systematic diflTerence in the calibrated results by taking the best 
20%, 50%, and 80% of all frames in terms of the fringe S/N 
(Tatulli et al. 2007 ). The diflTerence between the results obtained 
with the best 20% and 80% of frames is negligible, and there- 
fore, we took the best 80% of all frames for our final results. The 
errors of the resulting visibilities, DPs, and CPs were estimated 
in the manner described in Paper I. 

a Cen A (G2V) and a Cen B (KIV) were observed for the 
calibration of the interferometric data of Antares. a Cen is a 
triple system consisting of a Cen A, a Cen B, and Proxima 
Cen. However, we always had only one component (A or B) in 
the field of view of AMBER. We adopted angular diameters of 
8.314 ± 0.016 and 5.856 ± 0.027 mas derived by Kervella et al. 
(2003) for a Cen A and B, respectively. We also used a Cen A 
as a spectroscopic standard star to remove telluric lines from the 
observed spectra of Antares. However, the spectrum of a Cen A 
shows weak CO first overtone lines. In this case, as described 
in Ohnaka et al. (2012), the calibrated spectrum of the science 
target is derived as Ff = Fl^J(Fll\/Ff), where Ff^""^^^ 



^sci(cal) 
obs 



and 

denote the true and observed (i.e., including the atmo- 
spheric transmission and the detector's response) spectra of the 
science target (or calibrator), respectively. We used the high r es- 
olution solar spectrum presented by Wallace & Hinkle (119961 ) to 
estimate the true spectrum of the calibrator a Cen A, because 
this star has the same spectral type as the Sun. 

The wavelength calibration was done by using the telluric 
lines in the spectrum of a Cen A. As a template of the telluric 



Available at |http://www.jnimc.fr/data-processing_amber.htm| 
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Fig. 1. The uv coverage of our AMBER obser- 
vations of Antares in 2009 (a) and 2010 (b). 



lines, we convolved the atmospheric transmission spectra mea- 
sured at the Kitt Peak National Observator}0 to match the spec- 
tral resolutions of the data. The uncertainty in wavelength cali- 
bration is 1.7 X 10"^ yum (~2 km s"^). The wavelength scale was 
converted to the heliocentric frame and then to the laboratory 
frame using the heliocentric velocity of -3.5 km s"^ measured 
for Antares (Gontscharov 2006). Figure |2] shows the spectra of 
Antares observed in 2009 and 2010, together with the spectrum 
of Betelgeuse obtained in 2009 (Paper II). These spectra are 
binned to the same spectral resolution of 8000. The observed 
spectra of Antares are similar to that of Betelgeuse and show 
little time variation within an interval of one year. However, 
despite this absence of time variation in the observed spectra, 
our spectro-interferometric observations have detected signifi- 
cant time variation in the velocity field of the atmosphere, as 
presented in Sect. 13.31 



2.2. Image reconstruction 

We used the MiRA package ver.0.9.9EI(Thiebaut2008) to recon- 
struct an image from the interferometric measurements at each 
spectral channel. As described in Paper II, we first carried out 
the image reconstruction using simulated data. As described in 
Appendix ini we generated simulated stellar images (e.g., limb- 
darkened disk or stellar disk with inhomogeneities and/or an ex- 
tended component) and computed visibilities and CPs with re- 
alistic noise at the same uv points as in our AMBER observa- 
tions of Antares. With the true image known, we searched for the 
appropriate reconstruction parameters, such as the initial model 
and regularization scheme. We also applied the self-calibration 
technique described in Paper II, which restores the Eourier phase 
from the diflTerential phase measurements. 

3. Results 

3.1. Continuum data 

We selected the continuum spectral channels shortward of the 
CO band head at 2.295 yum, avoiding the CO lines as well as the 
weak lines. The average visibilities over the continuum spectral 
channels were computed in each data set. We adopted the sim- 
ple mean of the errors without reducing by V^cont, where A/cont 



^ http://www.eso.org/sci/facilities/paranal/instruments/isaac/tools/ 
spectra/atmos_S_K.fits 
^ http ://w w w-obs . univ-ly on 1 . fr/labo/per so/eric . thiebaut/mira. html 



is the number of the continuum spectral channels, because the 
errors are dominated by systematic errors, and they do not be- 
come smaller by averaging. 

Eigure[3]shows the continuum visibilities measured in 2009 
and 2010. The 2009 and 2010 data sample up to the fourth and 
fifth visibility lobes, respectively. This corresponds to angular 
resolutions of 9.8 and 7.2 mas for the 2009 and 2010 data, re- 
spectively, which mark the highest spatial resolution achieved 
for Antares at any wavelength. Uniform-disk fits to the 2009 and 
2010 data result in 36.97 ± 0.05 mas (reduced ;^^ = 4.6) and 
36.50 ± 0.05 mas (reduced ;^^ = 5.4), respectively. Eitting with 
a power-law-type limb-darkened disk (Hestroflfer et al. 1997 ) re- 
sults in a limb-darkened disk diameter of 37.38 ± 0.06 mas and 
a limb-darkening parameter of (8.7 ± 1.6) x 10"^ for the 2009 
data and 37.31 ± 0.09 mas and (1.5 ± 0.2) x 10"^ for the 2010 
data (the errors were estimated by the bootstrapping technique 
as described in Efron & Tibshirani 1 19931 1. 

The reduced ;^^ values for the fit to the 2009 and 2010 data 
are 3.9 and 3.3, respectively. These values are better than those 
with the uniform-disk fit but are still higher than 1. This is be- 
cause of the deviation from the limb-darkened disks in the data 
points at visibility nulls at spatial frequencies of -87 arcsec"^ 
(2009 data) and -115 arcsec"^ (2010 data), as well as those at 
95-100 arcsec"^ in the 2009 data. These last data points were 
obtained at position angles diflTering by 20-25°, which suggests 
the presence of inhomogeneities on a spatial scale of -10 mas in 
2009. However, as discussed in Sect. 13.31 the reconstructed im- 
ages in the continuum show only a very slight deviation from the 
limb-darkened disk: less than 2% (see Eig.[5}). The 2010 data in- 
clude uv points taken at position angles diflfering by roughly 90°, 
which are located at spatial frequencies of 130-138 arcsec"^ 
Still, Eig. [3] shows that the measured visibilities closely follow 
the limb-darkened disk. In summary, the overall deviation from 
the limb-darkened disks is small, indicating that the star shows 
only a weak signature of inhomogeneities in the continuum. 
Comparison between the 2009 and 2010 data suggests that a time 
variation with an interval of one year is only marginal. The same 
finding in Betelgeuse is reported in Paper II. Therefore, RSGs 
may show only a small deviation and time variation in the sur- 
face structure seen in the continuum. 

This observed time variation (or absence of it) is much 
smaller than the maximum visibility variation of ±40% pre- 
dicted by the current 3-D convection simulations in the third 
lobe for 2.2 yum, which approximately samples the continuum 
(see Eig. 18 of Chiavassa et al. 12009b . However, this does not 
necessarily mean that the observed time variation disagrees with 
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Fig. 2. Spectra of Antares observed in 2009 and 
2010, together with the spectrum of Betelgeuse 
from Paper II (the one obtained in 2009). 
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Fig. 3. The ^-band continuum visibilities of Antares observed in 2009 (red dots) and 2010 (blue dots). The fit with the power- 
law-type limb-darkened disk for the 2009 and 2010 data is shown by the red and blue solid lines, respectively. The limb-darkened 
disk has angular diameters of 37.38 and 37.31 mas and limb-darkening parameters of 0.087 and 0.15 for the 2009 and 2010 data, 
respectively. The insets show the enlarged views of the individual visibility lobes. 



the 3-D simulations. It is possible that we observed Antares at 
two epochs with small variations by chance, given that the 3-D 
simulations of Chiavassa et al. (2009) show that large convective 
cells have lifetimes of years. Furthermore, the standard devia- 
tion of the temporal variation in the visibility predicted by the 
3-D simulations is smaller than 40%, down to 10%, depending 
on the spatial frequency (see Figs. 11, 12, and 14 of Chiavassa 



et al. 120091) . It is also possible that significant inhomogeneities 
remained undetected due to the limited position angle coverage 
of our observations. Therefore, long-term monitoring AMBER 
observations with better position angle coverage are necessary 
to rigorously test the current 3-D convection simulations. 

The derived limb-darkened disk angular diameters are no- 
ticeably smaller than the 41.3+0.1 mas derived from lunar occul- 
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Fig. 4. H-R diagram with the theoretical evolutionary tracks for 
9, 12, 15, 20, and 25 M© stars for Z =Q.Q14 without rotation 
(solid lines) from Ekstrom et al. (120121 ) and the observationally 
derived position of Antares (filled diamond with error bars) . For 
the 15 and 9 Mq stars, the evolutionary tracks with rotation are 
plotted by the dashed lines. 



tation measurements of Richichi & Lisi (1990). However, they 
used a narrow-band filter at 2.43 yum with a FWHM of 0.035 yum. 
A number of CO first overtone lines and possibly also H2O lines 
from the MOLsphere are present in this wavelength range, which 
makes the star appear larger than in the continuum. Furthermore, 
Richichi & Lisi (1990) detected asymmetry in the brightness 
profile reconstructed from the lunar occultation data. This also 
indicates that their angular diameter is aflTected by the inhomo- 
geneous outer atmosphere, where the CO and H2O molecules 
form. 



3.2. Determination of basic stellar parameters 

By combining the above angular diameters in the continuum 
with the measured bolometric flux, we can determine the eflTec- 
tive temperature. We collected photometric data from the U band 
to 10 yum from the literature (Lee 1970 Ducati 2002 2MASS, 
Skrutskie et al. 2006). The photometric data were dereddened 
by the interstellar extinction, which was estimated as follows. 
The intrinsic B - V color of a red supergiant with the spec- 
tral type of Antares (M1.5Iab-b) is 1.70 (Elias et al. [T985]) . On 
the other hand, the observed B - V color of Antares is 1.84 
(Ducati 2002). Therefore, the color excess E{B - V) due to the 
interstellar extinction is 0.14, which translates into Ay = 0.43 
if Ay = 3AE(B - V) (Savage & Mathis ,1979.) is assumed. 



Using this interstellar extinction, we derived a dereddened bolo- 
metric flux of 8.36 X 10"^ W m"^. The irregular variability of 
Antares with an amplitude of ~ 1 magnitude in the visible leads 
to an uncertainty in the bolometric flux, because the photomet- 
ric data collected from the literature were not taken simultane- 
ously. However, the time variation in flux between the R and K 
bands is difl&cult to estimate from the data available in the liter- 
ature, because the data points are not sufl&cient. Therefore, we 
checked the time variation in flux in these bands for Betelgeuse 
and adopted it for Antares. We estimated the amplitude of the 
time variation in the R- and /-band flux of Betelgeuse to be 15% 
based on the photometric data from Ducati (2002) and Low et 
al. (1970). The light curves at 1.25, 2.2, 3.5, and 4.9 jum ob- 
tained over 3.6 years with COBE/DIRBIE suggest that the am- 
plitude of the flux variation is approximately 10% in these IR 
bands (Price et al. 2010). With these flux variations adopted for 
Antares, the uncertainty in the bolometric flux is estimated to be 
±1.07xlO-^Wm-2. 

The limb-darkened disk diameters derived from our AMBER 
continuum data and this bolometric flux result in an eflfective 
temperature of 3660 ± 120 K for both 2009 and 2010. This efl'ec- 
tive temperature agrees well with the eflTective temperature scale 
of Levesque et al. (2005), which gives 3710 K for the spectral 
type of Antares (Ml. 5). By combining the measured bolometric 
flux and the parallax of 5.89 ± 1.00 mas (van Leeuwen 2 0071) , 
we derive a luminosity of logL^/LQ = 4.88 ± 0.23. The linear 
radius estimated from the measured limb-darkened disk angu- 
lar diameters and the parallax is 680 Rq (3.2 AU). The eflTective 
temperature and luminosity we derived agree very well with Tqq 
= 370 7 + 77 K and logL^/L© = 4.99 ±0.15 derived by Pecaut 
et al. (I20T2I) . 

The observationally derived location of Antares on the H- 
R diagram is shown in Fig. [4] with the theoretical evolutionary 
tracks for 9, 12, 15, 20, and 25 Mq stars for the solar metal- 
licity Z=0.014 published in Ekstrom et al. (2012). The location 
of Antares agrees with the evolutionary tracks for a 15 Mq star. 
However, given the large uncertainty in the luminosity, the error 
in the initial stellar mass is estimated to be ±5 M©. This value 
agrees reasonably well with the 17.2 Mq derived by Pecaut et 
al. (2012), who estimated the mass using the evolutionary tracks 
with rotation of Ekstrom et al. (.2012) . The theoretical evolution- 
ary tracks of Ekstrom et al. (120121 ) shows that when a 15 Mq star 
reaches the observationally derived effective temperature and lu- 
minosity of Antares, the current mass is 13-14.3 M© (without 
rotation) and 11-14.3 Mq (with rotation). Combining the mass 
of 1 1-14.3 Mq and the linear radius of 680 Rq, we obtain a sur- 
face gravity of log g = -0.1. ..-0.2 (g is given in units of cm s"^ 
throughout the paper). 

The age of Antares is important for constraining the 
Upper Scorpius OB association, which is still controversial. 
Comparison of the position of Antares on the H-R diagram with 
the evolutionary tracks of Ekstrom et al. (120121 ) for a 15 Mq 
star with and without rotation suggests an age of 11-15 Myr. 
This agrees very well with the 12^^ Myr derived by Pecaut et 
al. (2012), which is not surprising given the good agreement for 
the eflTective temperature and luminosity between their results 
and ours. Pecaut et al. (120121) show that this age of Antares is 
consistent with the age of the Upper Scorpius OB association of 
9-13 Myr derived from stars with other spectral types. 

3.3. Reconstructed Images 

Figure [5] shows the reconstructed images from the 2009 data at 
four diflferent wavelengths (panels a-d) within the CO line ob- 
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Fig. 5. Reconstructed images of Antares in 2009 within one of the CO first overtone lines (2.3058 yum). a-d: The observed line 
profile with four representative wavelengths marked. Panels e-h show the images reconstructed at these wavelengths. The contours 
correspond to 1, 2, 4, 8, 16, 32, and 64% of the peak intensity of each image. Panel i shows the diflTerence between the reconstructed 
continuum image and the best-fit limb-darkened disk for the continuum visibilities, whose parameters are given in Sect. 13.11 Panels 
j-1 show the continuum- subtracted images, in which the continuum image (panel e) is subtracted from the images in the panels f-h. 
The limb of the star as defined in Sect. l3.3l is marked by the red solid lines. 



served at 2.3058 yum0. All reconstructed images appear strongly 
elongated due to the elongated beam (9.8 x 63.2 mas), which re- 
sults from the lack of uv points roughly in the North-South direc- 
tion (see Fig. [T]). The reconstructed continuum image (Fig. [5^) 
shows only a slight deviation from the limb-darkened disk that 
represents the best fit to the continuum visibilities. This is 
clearly seen in Fig. [5j, which shows the residual continuum 
image obtained after subtracting the limb-darkened disk with 
the parameters derived in Sect. 13.1 [ convolved with the interfer- 
ometer's beam (the determination of the beam is described in 
Appendix|B]). The residuals are at most 2% of the peak intensity 
of the reconstructed continuum image. 

The image in the blue wing (Fig.fSf) appears very similar to 
that in the continuum except for the very weak extended com- 
ponent on either side of the star. However, the intensity of this 
extended component is comparable to the noise level of the im- 
age reconstruction (~1%), which is estimated from the intensity 



^ A movie of the data cube containing all wavelengths from 2.28 and 
2.31 jum is available at 



http://www.mpifr-bonn.mpg.de/staff/kohnakal 



of the background features sufl&ciently away from the star in the 
reconstructed images. Therefore, it is not clear whether or not 
this extended component in the blue wing image is real. On the 
other hand, the images in the line center and red wing (Figs. |5^ 
and[5]i) show a much more prominent, asymmetrically extended 
component. 

In order to clearly show the diflTerence among the images, we 
normalized each image with its peak intensity and subtracted the 
normalized continuum image from the images within the line. 
These continuum- subtracted images are shown in the bottom 
row of Fig. [5] Also plotted is the limb of the star (red solid line), 
which we define by the contour at 3% of the peak intensity of the 
continuum image. We adopted this contour because the CO lines 
appear in emission outside the limb defined in this manner as de- 
scribed below (Fig. [6]). Note that we do not artificially adjust the 
position of the images (either in the lines or in the continuum) at 
all, because the relative astrometry among the images at diflTer- 
ent spectral channels is conserved thanks to the self-calibration 
technique (see Paper II for details). The continuum- subtracted 
images in the line center and red wing clearly show an extended 
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component that is more prominent on the western side of the star, 
while the continuum- subtracted image in the blue wing shows 
almost no extended component. The resolved extended compo- 
nent is present off the limb of the star (i.e., outside the stellar 
disk) particularly on the western side of the line center image 
and also-albeit weaker-in the red wing image. This represents 
the second imaging of the MOLsphere in RSGs within individ- 
ual CO lines after Betelgeuse reported in Paper II. 

While we normalized each of the above images with its peak 
intensity, it is also possible to normalize each image so that the 
intensity integrated over the entire image is equal to the observed 
flux at that spectral channel. This allows us to extract the spa- 
tially resolved spectrum at each position over the stellar image. 
The spatially resolved 2009 spectra at three representative posi- 
tions along the direction of the beam's minor axis are shown in 
Fig. [6j which reveals the CO emission lines ofl' the limb on the 
western side (position AJj. These emission lines are redshifted 
with respect to the absorption lines in the spatially unresolved 
observed spectrum, reflecting the appearance of the extended 
component only in the line center and red wing. On the east- 
em side, the line profiles are characterized by weak blueshifted 
emission and redshifted absorption. This is explained as follows. 
Because of the finite beam size, we see not only the emission oflT 
the limb but also the absorption originating from inside the limb. 
The weak blueshifted emission, which is seen as the weak ex- 
tended feature in the reconstructed blue wing image (Fig. |5f), 
partially fills in the absorption. This results in what appears as 
redshifted absorption with weak blueshifted emission. Ohnaka 
(12013b presents similar spatially resolved CO first overtone line 
spectra for Betelgeuse and reveals the CO line emission off the 
limb, but for the 1-D image. Our 2-D image reconstruction for 
Antares demonstrates that it is now feasible to obtain the spa- 
tially resolved spectrum at each position over the surface of a 
star as well as oflf the limb of the star, as is routinely done in 
solar physics. 

Figure |7] shows the reconstructed images from the 2010 data 
for the same CO line as shown in Fig. [S] Because the position 
angle coverage of the 2010 observations is better than the 2009 
observations, the beam of the 2010 images is less elongated than 
that of the 2009 images. However, for a direct comparison of the 
images taken at two epochs, the 2010 images are convolved with 
the same beam as used for the images in 2009. The reconstructed 
image in the continuum shows little deviation from the best- fit 
limb-darkened disk derived from the continuum visibilities, as in 
the 2009 data. The residual after subtracting the limb-darkened 
disk with the parameters derived in Sect. 13 . 1 1 is at most 1 % of 
the peak intensity of the reconstructed continuum image. The 
continuum- subtracted images shown in the bottom row of Fig. [7] 
reveal that the images in the blue wing and line center show a 
prominent extended component. On the other hand, the image 
in the red wing shows little trace of this (the intensity of the 
extended component in the red wing image is ~2%, which is 
only slightly higher than the image reconstruction noise of 1%). 
This means that the appearance of the star in the blue and red 
wing has swapped within one year (compare the images in the 
bottom row of Figs. [5] and (T]). 

Figure[8]shows the spatially resolved spectra at three selected 
positions over the stellar image. Obviously, the blueshifted emis- 
sion appears oflf the limb (position A), although the strength of 
the emission is weaker than the redshifted emission seen in 2009. 



^ A movie of the spatially resolved spectra across the stellar image is 

availableat 

|http://www.mpifr-bonn.mpg.de/staff/kohnaka| 
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Fig. 6. Spatially resolved spectra obtained in 2009 at three rep- 
resentative positions over the stellar image are plotted by the red 
solid lines. The representative positions are marked in the top 
panel, which shows the line center image at 2.30580 yum (same 
as Fig. [5^). The black lines represent the spatially unresolved 
spectrum. The spatially resolved spectra are scaled for compari- 
son with the spatially unresolved spectrum. The black solid line 
in the top panel marks the limb of the star. 



This is simply because the extended component, MOLsphere, is 
weaker in 2010 than in 2009. This can also be seen in the visi- 
bilities observed in the first lobe, as shown in Fig. [9] While the 
observed spatially unresolved spectra show little time variation, 
the visibilities in the CO lines observed in 2010 are noticeably 
higher than observed in 2009. This means that the star appears 
less extended in the CO lines in 2010 than in 2009. The figure 
also shows the time variation in the visibilities within the CO 
lines, which corresponds to the diflTerent appearance of the star 
within the CO lines as shown in Figs. [5] and [TJ 

The time variation in the reconstructed images within the 
CO lines suggests a significant change in the dynamics of the 
atmosphere. Such a time variation in the atmospheric dynamics 
has already been detected in Betelgeuse (Paper II). Therefore, 
our AMBER observations of Antares implies that the inhomoge- 
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Fig. 7. Reconstructed images of Antares in 2010 shown in the same manner as Fig. [5] 



neous, highly temporally variable nature of the outer atmosphere 
might be common among RSGs. 



4. Modeling of the velocity field 

The diflTerent appearance of the star within the individual CO 
lines was also detected in Betelgeuse, and it was interpreted as 
inhomogeneous gas motions in the photosphere and MOLsphere 
(Papers I and II). To characterize the velocity field using the 
above reconstructed images of Antares, we applied the same 
two-layer model as used for Betelgeuse (see Paper I for details 
of the model). The star is represented with the blackbody of 
3700 K, which is the eff'ective temperature of Antares derived 
from our continuum visibilities (Sect. [T2] ). We represent the CO 
gas in the photosphere and MOLsphere with an inner and outer 
layer, respectively. The reason we adopted this two-layer model 
is that we cannot reproduce the observed CO line spectrum and 
the CO line images simultaneously with a single layer model. 
For example, single layer models with CO column densities of 
10^^-10^^ cm"^ and a radius of ~ 1.5 7?^ can roughly explain 
the observed extended emission in the CO lines. However, these 
models predict the CO absorption spectrum to be too weak com- 
pared to the observations or to appear even in emission, because 
of the strong extended emission from the layer. This suggests 
that the inner layer representing the geometrically compact pho- 



tosphere is necessary to explain the images and spectrum simul- 
taneously. 

Each layer is characterized by its radius, CO column density, 
and temperature. As in Papers I and II, the layers are assumed to 
be geometrically thin, and the thickness of the layers are fixed to 
0.02 Ri,. Because this choice is arbitrary, we computed models 
with larger geometrical thicknesses (up to 0.05 and 0.2 R^, for 
the inner and outer layer, respectively) and confirmed that the 
thickness of the layers does not aff'ect the results of the model- 
ing. Changing the inner layer radius aff'ects the intensity distribu- 
tion only slightly. The spatial resolution of the current AMBER 
data is sufficient for constraining the size the MOLsphere (i.e., 
outer radius) but insufficient for constraining the inner radius. 
We approximate an inhomogeneous velocity field with a patch 
(or clump) of CO gas that is moving at a diff'erent velocity from 
the gas in the remaining region. 

For the inner CO layer, we adopted the same parameters as 
for Betelgeuse, because the eff'ective temperature and surface 
gravity of Antares are similar to those of Betelgeuse (Tqq = 
3690 K from Paper II and log g = -0.3 from the mass and radius 
estimated by Harper et al. 120081 ). We set the radius, temperature, 
and CO column density of the inner layer to be L05 R^,, 2250 K, 
and 5 x 10^^ cm"^, respectively, as in Papers I and II. The ra- 
dius, temperature, and CO column density of the outer layer, as 
well as the velocity field, were treated as free parameters, and 
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Fig. 8. Spatially resolved spectra obtained in 2010 shown in the 
same manner as Fig. (6] 



we assumed the same velocity field for the inner and outer lay- 
ers. Because the elongated beam of our AMBER observations 
makes it difi&cult to know the actual number and shape of stellar 
spots and the spatial resolution in the North-South direction is 
low, we assumed only a single patch of CO gas to keep the num- 
ber of free parameters as small as possible. The model images 
were spectrally convolved with the spectral resolution of the ob- 
served data (A/AA = 8000) and then convolved with the same 
beam as for the reconstructed images. 

Figures [TOf-h show the best-fit continuum- subtracted model 
images derived for the 2009 observations within the same CO 
line as shown in Fig. [5] A comparison between the observed line 
profile and that predicted by the model is also shown. This model 
is characterized by a single large patch (or clump) of CO gas in- 
falling with 2.5 km s"^ which is marked with the dashed circle 
in the red wing image in Fig. [TOb (see also the schematic view 
in Fig. O. This CO gas patch or clump produces slightly red- 
shifted absorption, which explains why the intensity within the 
patch in the red wing image (Fig. [TOb ) is lower than in the re- 
maining region. The CO gas outside the patch is moving outward 
much faster at 25 km s"\ as marked in the line center image in 
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Fig. 9. Comparison of the visibilities observed in 2009 (data set 
#26, black solid line) and 2010 (data set #107, red solid line). 
Both visibilities were taken at almost the same projected base- 
line length (13.9 m) and position angle (84°-85°). The spatially 
unresolved spectra observed in 2009 and 2010 are also plotted 
by the green and blue solid lines. 



Fig.fTOb (see also the schematic view in Fig.[T2|. The outer CO 
layer extends to 1.3 Ri, with a temperature of 2000 K and a CO 
column density of 1 x 10^^ cm"^. The extended component ob- 
served in the line center and red wing, as well as the absence of 
it in the blue wing, is reasonably reproduced. From the observed 
images, the velocity of the CO gas moving inward and outward 
is constrained to be 0-5 km s"^ and 20-30 km s"\ respectively. 
However, the model predicts the extended component to be too 
strong, and the observed asymmetry of the extended component 
(i.e., it extends only to the western side) is not reproduced. The 
model flux is also too narrow compared to the observed line pro- 
file. 

A summary of our modeling is given in Table \T\ We also 
computed the reduced ;^^ value for each model in the continuum, 
blue wing, line center, and red wing from the fit to the observed 
visibilities and closure phases. The aforementioned diflferences 
between the model and the observed images are reflected in the 
reduced x^ values that are much higher than 1 . This is presum- 
ably due to the simplifications adopted in our model. In partic- 
ular, we assumed the (radial) column density of CO and tem- 
perature in each layer to be constant regardless of the direction 
with respect to the stellar center and attempted to explain the ob- 
served results only by introducing a simplified velocity field. We 
also assumed a single patch to keep the number of free parame- 
ters as small as possible. Therefore, inhomogeneities in the CO 
density and temperature, as well as multiple patches, may recon- 
cile the disagreement between the observed data and the model. 
However, our goal is to obtain an approximate picture of the ve- 
locity field, and therefore, we refrain from introducing more free 
parameters, which cannot be well constrained from the present 
data. 

The best-fit continuum- subtracted model images for the 
2010 data are shown in Fig. \TT\ In this model, a large patch 
of CO gas is moving outward only with 2.5 km s"^ while the 
CO gas in the remaining region is infalling much faster with 
25 km s"^ (see also the schematic view in Fig. [12]). The patch 
is so large that it dominates the surface of the star facing to- 
ward the observer. The radius and temperature of the outer CO 
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Fig. 12. Schematic view of the velocity fields in Antares derived 
from our AMBER observations in 2009 and 2010. Only the outer 
CO layer is drawn for the sake of clarity of the figure, but the 
velocity field of the inner CO layer was assumed to be the same 
as in the outer layer The images are not to scale. 



layer are the same as those for the 2009 data, but the CO col- 
umn density (5 x 10^^ cm"^) is smaller than that derived from 
the 2009 data. The extended component observed in the line 
center and blue wing is fairly reproduced, and the model im- 
age in the red wing shows little trace of the extended component 
as observed. The model predicts the extended component to be 
too strong compared to the observed image, and the line profile 
predicted by the model is also too shallow. However, as men- 
tioned above, this disagreement is likely due to the simplifica- 
tions adopted in our model. The smaller CO column density of 
the outer CO layer found for the 2010 data reflects the weaker 
extended component in 2010 compared to 2009. However, we 
note that this weaker extended component can also be explained 
by a model with a smaller radius of 1.2 R^, and a CO column 
density of 1 x 10^^ cm"^ for the outer CO layer We also found 
that the models with the up welling CO gas with 0-5 km s"^ and 
downdrafting gas with 20-30 km s"^ can explain the observed 
data. The uncertainties in the radius, temperature, and CO col- 
umn density of the outer CO layer are estimated to be ±0.1 7?^, 
±200 K, and a factor of 2, respectively. 

Comparison between the velocity field obtained from the 
2009 and 2010 data reveals that the overall gas motions in the 
outer atmosphere changed the direction within one year: the 
strong up welling motions were dominant in 2009, while the ve- 
locity field in 2010 was dominated by strong downdrafts. 



5. Discussion 

The characteristics of the velocity field found in Antares are sim- 
ilar to those reported for Betelgeuse in Papers I and II. In both 
stars, we found CO gas motions with velocity amplitudes of up 
to 20-30 km s"^ on a spatial scale larger than the stellar angular 
radius. Significant time variations in the velocity field have been 
detected over an interval of one year The radius, temperature, 
and density of the outer CO layer derived for Antares are also 
similar to those in Betelgeuse. Therefore, these inhomogeneous 
velocity fields in the MOLsphere extending to 1.2-1.5 7?^ might 
be common among RSGs. 

A schematic view of the derived velocity fields in 2009 and 
2010 is given in Fig. [121 It is not yet known whether all out- 
wardly moving gas clumps fall back after some time or they are 
further accelerated outward. If we assume that the inhomoge- 
neous CO gas motions in Antares are ballistic, the CO gas mov- 
ing outward with 20-30 km s"^ seen at 1.3-1.4 R^, in April 2009 
reaches a maximum height of 1.45-1.60 7?* in 2-3 months and 
falls back to the original height in another 2-3 months. Likewise, 



'-30 km/s if we assume that the CO gas infalling fast at 20-30 km s seen 



in 2010 is falling back from some height with an initial veloc- 
ity of zero, it is estimated to have traveled from 1.45-1.60 Ri, 
over 2-3 months. On the other hand, the CO gas moving inward 
at only 0-5 km s"^ (2009) or outward (2010) may represent a 
clump that has traveled from a deeper layer and is near the max- 
imum height, at which the velocity is zero. In this case, the clump 
is estimated to have traveled from -1.2 7?* with an initial veloc- 
ity of 20-30 km s"^ over approximately two months. This means 
that multi-epoch AMBER observations over 2-6 months would 
be useful for probing whether the outwardly moving gas clump 
falls back or it is accelerated outward further 

The physical mechanism responsible for the formation of the 
MOLsphere and the inhomogeneous gas motions in this region is 
still unclear The upwelling and downdrafting motions are rem- 
iniscent of convection, which may levitate the gas to form the 
MOLsphere. For example, Lim et al. (119981 ) suggest that con- 
vection could lift the photospheric material into the extended at- 
mosphere. If convection alone is responsible for the formation of 
the MOLsphere, 3-D convection simulations should predict the 
density consistent with the observationally derived values at the 
location of the MOLsphere (1.2-1.4 R^). 

Given that the current 3-D convection simulations cover up 
to 8-9 Ri,, it is possible to carry out this comparison. One of 
the 3-D models of Chiavassa et al. (12011b , which is closest to 
Antares or Betelgeuse, shows that the density at 1.2 7?^ ranges 
between 10"^^ and 10"^^ g cm"^ (see their Fig. 4). On the other 
hand, we can estimate the gas density in the MOLsphere as fol- 
lows. While we adopted the geometrical thickness of 0.02 7?^ 
for the outer CO layer representing the MOLsphere, the in- 
ner radius of the MOLsphere is not observationally constrained. 
Therefore, we estimated a lower limit on the CO density by di- 
viding the CO column density with the geometrical thickness 
defined by the outer radius of the MOLsphere (i.e., geometrical 
thickness = 0.2-0.4 R^,). If we assume chemical equilibrium in 
the MOLsphere, we can estimate the gas density from the CO 
density and the temperature of the MOLsphere (2000 K). The 
derived gas density is ~ 10"^^ cm"^, which is higher than the 
values predicted by the 3-D convection model by six to eleven 
orders of magnitude. The analysis of line profiles in Antares by 
Gray & Pugh (2012) also implies that convective cells pene- 
trate only the lower portion of the photosphere, although it is not 
straightforward to convert the radial velocities derived from spa- 
tially unresolved spectra to the actual velocity field. Therefore, 
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Fig. 10. Best-fit stellar patch model with an inhomogeneous velocity field for Antares observed in 2009. Panels a-d show model 
images in the continuum, blue wing, line center, and red wing within the same CO line profile as shown in Fig. [5] The wavelengths 
of these images are marked with black, blue, green, and red dots in Panel e, where the observed and model spectra are plotted by the 
filled diamonds and the black solid line, respectively. Panels f-h show the best-fit continuum- subtracted model images convolved 
with the same beam as the observed images. Therefore, these model images can be compared to the images shown in Figs.[5]-[5j. 
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Fig. 11. Best-fit stellar patch model with an inhomogeneous velocity field for Antares observed in 2010, shown in the same manner 
as Fig. [TOl The continuum- subtracted images from the best-fit model, which are shown in Panels f-h, are compared to the images 
shown in Figs.[7]-[7j. 
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Table 1. Parameters of the MOLsphere model with an inhomogeneous velocity field for Antares. The definition of the patch size 
and position is depicted in Fig. 4 of Paper I. 



Parameter 


Value (2009) 


Value (2010) 


Inner layer 


CO column density (cm""^) 
Temperature (K) 
Radius (R^) 


5 X 10"" (fixed) 
2250 (fixed) 
1.05 (fixed) 


5 X 10"" (fixed) 
2250 (fixed) 
1.05 (fixed) 


Outer layer 


CO column density (cm""^) 
Temperature (K) 
Radius (R^) 


1 X 10'^ (+0.3 dex) 
2000 + 200 
1.3 + 0.1 


5 X 10^^ (+0.3 dex) 
2000 + 200 
1.3 + 0.1 


Velocity field 


Velocity within the patch (km s~^) 
Velocity outside the patch (km s~^) 
Stellar patch size (0)(°) 
Stellar patch position (6, cp) (°) 


2.5 + 2.5 (inward) 
25 + 5 (outward) 
60+10 
(0-30, 0-360) 


2.5 + 2.5 (outward) 
25 + 5 (inward) 
80+10 
(0-30, 0-360) 


Reduced ;^'^ 


Continuum 
Blue wing 
Line center 
Red wing 


10.6 
32.6 
21.0 
29.8 


18.1 
31.4 
18.9 
20.5 



at the moment, convection alone is unlikely to explain the for- 
mation of the MOLsphere and the inhomogeneous gas motions. 
However, it remains to be tested whether the combination of con- 
vection and some other mechanism (e.g., radiation pressure on 
molecules and atoms or scattering on dust grains) can explain 
the formation of the MOLsphere and the mass loss in RSGs. 

Another possible mechanism is the propagation of Alfven 
waves. Magnetic fields have been detected in Betelgeuse 
(Auriere et al. 2010). While Grunhut et al. (2010) detected no 
magnetic fields in Antares, this can be due to the insufl&cient 
S/N as they themselves note. Therefore, future observations with 
a better S/N may reveal magnetic fields in more RSGs including 
Antares, and the magnetohydrodynamical (MHD) wave propa- 
gation remains a possible process to explain the formation of 
the MOLsphere and the inhomogeneous gas motions. However, 
the currently available MHD model by Airapetian et al. (120001) 
deals with the stellar winds from the chromosphere and does 
not include the MOLsp here. T he recent magneto-rotator model 
by Thirumalai & Heyl (2012) represents an interesting alterna- 
tive for the RSG mass loss, although their model is still station- 
ary and does not include the chromosphere. The inclusion of 
a proper treatment of the time-dependent, multicomponent na- 
ture (chromosphere and MOLsphere) would be crucial for test- 
ing whether or not the MHD process with or without rotation is 
responsible for the formation of the MOLsphere and the inho- 
mogeneous gas motions. 

While the presence of the MOLsphere is confirmed by the 
IR interferometric imaging of Antares and Betelgeuse in the in- 
dividual CO lines, there are still problems with the MOLsphere, 
as mentioned in Ohnaka et al. (120 121) . The additional TiO ab- 
sorption originating in the MOLsphere makes the TiO bands 
too strong compared to the observed spectra of RSGs (Hron et 
al. 2010). As Ohnaka et al. ("2012) suggest, the scattering due 
to TiO molecules may reconcile the above problem, because 
the line formation by scattering can be important for electronic 
transitions such as the TiO bands in the visible. Another prob- 
lem is that the current MOLsphere models cannot explain the 
high resolution 12 yum spectra of H2O lines (Ryde et al. 2006). 
Recently, Lambert et al. (2013) suggest the importance of NLTE 
eflTects in the formation of these H2O lines. Therefore, model- 
ing with the scattering of molecules and NLTE eflfects taken into 



account is necessary for a comprehensive understanding of the 
MOLsphere. 

6. Concluding remarks 

We have presented high spectral resolution aperture- synthesis 
imaging of Antares at two epochs with VLTI/ AMBER. With the 
stellar surface spatially resolved with AMBER's high spatial res- 
olution, amber's high spectral resolution has enabled us to de- 
rive the kinematical information at different positions over the 
surface of the star. The reconstructed images obtained at the first 
epoch (April 2009) reveal an asymmetrically extended compo- 
nent (MOLsphere) in the red wing and line center of the individ- 
ual lines, while the image in the blue wing shows little trace of 
the extended component. At the second epoch (April 2010), the 
extended component appears in the blue wing and line center, 
and the red wing image does not exhibit the extended compo- 
nent. We have also found that the extended component became 
weaker in 2010. These results can be interpreted as a significant 
time variation in the atmospheric dynamics. 

Our modeling suggests that the velocity field in the 
MOLsphere in 2009 is characterized by a large patch or clump 
of CO gas infalling at only 0-5 km s"^ with the CO gas in the re- 
maining region moving outwardly much faster at 20-30 km s"^ . 
The data observed in 2010 are explained by a model in which 
a large patch of CO gas moving outwardly at 0-5 km s"^ and 
the CO gas in the remaining region is infalling faster at 20- 
30 kms"^. Our modeling also shows that the MOLsphere is 
extended to 1.2-1.4 R^, with column densities of 5 x 10^^- 
1x10^^ cm"^ and a temperature of -2000 K. These properties of 
the MOLsphere of Antares, as well as the velocity field, are very 
similar to those derived for Betelgeuse. Therefore, the inhomo- 
geneous gas motions in the MOLsphere extending to 1.2-1.5 7?^ 
might be common among RSGs. 

The origin of the detected gas motions is not clear yet. To 
understand the physical mechanism responsible for the dynam- 
ics and formation of the MOLsphere, it is necessary to spatially 
resolve the dynamics at diff'erent heights from the deep photo- 
sphere to the outer atmosphere. Similar AMBER observations 
for diff'erent atomic and molecular lines would make this pos- 
sible. For example, weak atomic lines such as Ti, Na, and Ca 
would be useful for probing the dynamics in the deep photo- 
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spheric layers. Moreover, it is essential to obtain images with 
better uv coverage. This is now feasible with more telescope con- 
figurations appropriate for the aperture- synthesis imaging avail- 
able at the VLTI. 
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Appendix A: Summary of AMBER observations 

Our AMBER observations of Antares are summarized in 
Tables|A3]and|A2l 



Appendix B: Image reconstruction of simulated 
data 

We generated interferometric data (visibility amplitudes and clo- 
sure phases) from simulated stellar images by sampling at the 
same uv points as our AMBER observations. The image recon- 
struction from these simulated interferometric data allows us to 
find the appropriate parameters (such as the initial model, prior, 
and regularization scheme) that can restore the original images 
correctly. 

For our data on Antares, we tested the image reconstruc- 
tion for a limb-darkened disk with the angular diameter and the 
limb-darkening parameter derived for the 2009 data and a uni- 
form disk with surface inhomogeneities and an extended compo- 
nent. The simulated interferometric data were generated using a 
program developed by one of the authors (K.-H. Hofmann) as 
follows. For each telescope triplet, we simulated two-telescope 
interferograms corresponding to the three baseline vectors for 
a given object (e.g., limb-darkened stellar disk or spotted star). 
The simulated interferograms were degraded by the atmospheric 
piston, photon noise, sky background, and detector noise. From 
-1000 simulated two-telescope interferograms on each base- 
line, the average power spectrum and bispectrum were calcu- 
lated. The subtraction of the noise bias terms and calibration 
with an unresolved calibrator star yielded the calibrated visibili- 
ties and closure phases for the simulated object. The amount of 
noise in the simulated interferograms was chosen to obtain ap- 
proximately the same errors as in the AMBER measurements of 
Antares. 

The reconstructed images were convolved with the clean 
beam, which was determined in the following manner. The spa- 
tial resolution in the direction of the uv points with the longest 
baseline length (^max) is given by /l/^max = 9.8 mas. However, 
the sparse uv coverage shown in Fig. [T] leads to a strongly elon- 
gated beam. We fitted the central peak of the dirty beam with a 
2-D Gaussian and derived the ratio between the major and minor 
axes and the position angle. The beam defined in this manner is 
narrower than that given by /l/^max, because of the lack of data at 
very short baselines. Therefore, we scaled the major and minor 
axes of this Gaussian beam so that the beam size along the mi- 
nor axis matches the resolution given by /l/^max- The final clean 
beam for the uv coverage obtained in 2009 is a 2-D Gaussian 
with 9.8 X 63.2 mas (FWHM) with the major axis at a position 
angle of -12°. 

The results of the image reconstruction of two sets of sim- 
ulated data are shown in Figs. IB. II and IB. 21 The reconstructed 
images agree reasonably well with the original image after con- 
volving with the beam. These tests show that the uniform-disk 
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Table A.l. Log of AMBER observations of Antares and the calibrator a Cen A in April 2009 with the EO-GO-HO (16-32-48m) 
baseline configuration. The seeing (s) and the coherence time (tq) are in the visible. 
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(m) 



PA 

n 



s 

n 



To 

(ms) 



NDIT 



(UTC) 



^P 
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PA 

n 



s 

n 
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Antares 
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2009 Apr 25 


1 


05:14:38 
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58/ 58/ 58 


0.65 


4.2 


500 


36 
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56/ 56/ 56 


0.52 


5.0 


200 


2 


05:16:57 
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58/ 58/ 58 


0.79 


3.6 


500 


37 
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0.53 


4.9 
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3 
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Of Cen A 



2009 Apr 24 



2009 Apr 25 



CI 
C2 
C3 
C4 
C5 



05:35:10 
06:48:04 
07:31:50 
08:10:33 
08:51:55 



30.98/15.50/46.48 
29.38/14.70/44.07 
28.10/14.06/42.16 
26.84/13.43/40.27 
25.46/12.74/38.20 



81/81/81 

96/ 96/ 96 

107/107/107 

116/116/116 

127/127/127 



0.74 
0.77 
0.40 
0.43 
0.51 



3.8 
3.3 
6.4 
5.8 
5.0 



2500 
2500 
2500 
2500 
2500 



C6 05:30:35 

C7 06:23:17 

C8 07:56:55 

C9 08:41:46 

CIO 09:45:10 



30.99/15.51/46.49 
29.91/14.97/44.87 
27.16/13.59/40.75 
25.66/12.84/38.50 

1 135.51 



80/ 80/ 80 

92/ 92/ 92 
114/114/114 

125/125/125 
-/-/144 



0.65 
0.60 
0.43 
0.35 



4.0 2500 

4.2 2500 

7.3 2500 
9.7 2500 



0.61 6.8 2500 



initial model with angular diameters of 36-37 mas and the same 
prior as used in Paper II are appropriate for the image recon- 
struction from our AMBER data of Antares. The prior is given 
by 



Appendix C: Fit to tlie interferometric data 

Figures lCTI and IC2l show the fit to the measured interferometric 
observables for the image reconstruction in the CO line shown 
in Figs. [5] and [71 respectively. 



Pr{r) 



1 



^(r-rp)/ep _|_ ^ ' 



where r is the radial coordinate in mas, and r^ and s^ define 
the size and the smoothness of the edge, respectively {s^ -^ 
corresponds to a uniform disk). The appropriate values for r^ 
and £^ were found to be 12 mas and 2.4 mas, respectively. We 
adopted the maximum entropy regularization as in Paper II. 
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Fig. B.l. Image reconstruction from the simulated data for a limb-darkened disk with the parameters derived from the 2009 data. 
a: Original image, b: Original image after convolving with the Gaussian beam with 9.8 x 63.2 mas. c: Reconstructed image after 
convolving with the beam, d: The intensity profiles of the original image are plotted by the red dotted line (along the minor axis) 
and the black dotted line (along the major axis). The intensity profiles of the reconstructed images are plotted by the red solid line 
(along the minor axis) and black solid line (along the major axis). The abscissa for the intensity profiles along the major axis is 
shown at the top of the panel, e: Comparison between the visibilities generated from the original image (filled dots) and those from 
the reconstructed image (triangles), f: Comparison between the CPs generated from the original image (filled dots) and those from 
the reconstructed image (triangles). 
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Fig. B.2. Image reconstruction from the simulated data for a uniform disk with surface inhomogeneities and an asymmetric extended 
component, shown in the same manner as Fig. IB. II 
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Table A.2. Log of AMBER observations of Antares and the calibrators a Cen A and a Cen B in April 2010. The data sets #66-99 
and C#l 1-C#18 were taken with the EO-GO-HO (16-32-48m) telescope configuration, while the data sets #100-109 and C#19 were 
taken with the EO-GO-Il (16-57-69m) configuration. The seeing (s) and the coherence time (tq) are in the visible. 



(UTC) 



(m) 



PA 

n 



s 

n 



To 

(ms) 



NDIT 



(UTC) 



(m) 



PA 

n 



s 

n 



To NDIT 
(ms) 



Antares 



2010 Apr 10 


2010 Apr 11 


66 


06:19:19 


45.55/15.19/30.36 


59/ 59/ 59 


1.02 


47.1 


500 


89 


05:47:04 


43.95/14.66/29.30 


55/55/55 


0.78 


30.5 


500 


67 


06:21:30 


45.66/15.23/30.44 


59/ 59/ 59 


0.97 


31.3 


500 


90 


05:49:16 


44.09/14.70/29.39 


55/55/55 


0.74 


23.9 


500 


68 


06:23:41 


45.77/15.26/30.51 


60/ 60/ 60 


0.99 


28.1 


500 


91 


05:51:27 


44.22/14.75/29.47 


55/55/55 


0.75 


23.1 


500 


69 


06:25:51 


45.88/15.30/30.58 


60/ 60/ 60 


1.11 


25.9 


500 


92 


05:53:37 


44.35/14.79/29.56 


56/56/56 


0.67 


30.7 


500 


70 


06:28:02 


45.98/15.33/30.65 


60/ 60/ 60 


1.07 


29.7 


500 


93 


05:55:48 


44.48/14.83/29.65 


56/56/56 


0.69 


27.5 


500 


71 


07:16:22 


47.64/15.89/31.76 


67/ 67/ 67 


0.83 


21.8 


500 


94 


06:34:20 


46.44/15.49/30.95 


62/62/62 


0.72 


17.1 


500 


72 


07:18:32 


47.69/15.90/31.78 


67/ 67/ 67 


0.81 


21.8 


500 


95 


06:36:31 


46.53/15.52/31.01 


62/62/62 


0.69 


17.8 


500 


73 


07:20:43 


47.73/15.92/31.81 


67/ 67/ 67 


0.78 


20.5 


500 


96 


06:38:41 


46.61/15.54/31.07 


62/62/62 


0.68 


18.2 


500 


74 


07:22:53 


47.77/15.93/31.84 


68/ 68/ 68 


0.84 


17.6 


500 


97 


06:40:52 


46.70/15.57/31.13 


63/63/63 


0.66 


19.5 


500 


75 


07:25:03 


47.80/15.94/31.86 


68/ 68/ 68 


0.85 


18.0 


500 


98 


06:43:03 


46.78/15.60/31.18 


63/63/63 


0.67 


20.3 


500 


76 


08:44:42 


46.97/15.66/31.30 


77/ 77/ 77 


0.68 


15.4 


500 


99 


07:20:36 


47.79/15.94/31.86 


68/68/68 


0.70 


34.6 


500 


77 


08:46:53 


46.88/15.64/31.25 


77/ 77/ 77 


0.66 


15.7 


500 


100 


09:21:33 


65.11/54.08/14.96 


118/128/81 


0.75 


20.6 


500 


78 


08:49:03 


46.80/15.61/31.19 


77/ 77/ 77 


0.69 


15.0 


500 


101 


09:23:43 


64.94/53.97/14.91 


119/128/81 


0.74 


21.7 


500 


79 


08:51:14 


46.71/15.58/31.13 


77/ 77/ 77 


0.71 


14.8 


500 


102 


09:25:54 


64.77/53.85/14.86 


119/129/81 


0.76 


23.1 


500 


80 


08:53:25 


46.61/15.55/31.07 


77/ 77/ 77 


0.69 


14.8 


500 


103 


09:28:05 


64.60/53.74/14.81 


119/129/81 


0.78 


23.4 


500 


81 


09:36:33 


44.08/14.70/29.38 


82/ 82/ 82 


0.70 


15.2 


500 


104 


09:30:15 


64.42/53.62/14.76 


120/129/81 


0.78 


23.4 


500 


82 


09:38:45 


43.92/14.65/29.27 


82/ 82/ 82 


0.85 


12.8 


500 


105 


09:56:49 


62.04/52.07/14.04 


124/134/84 


0.71 


30.4 


500 


83 


09:40:55 


43.75/14.59/29.16 


82/ 82/ 82 


0.88 


12.8 


500 


106 


09:58:59 


61.83/51.93/13.98 


124/134/84 


0.68 


28.1 


500 


84 


09:43:06 


43.58/14.54/29.05 


82/ 82/ 82 


0.85 


12.7 


500 


107 


10:01:09 


61.61/51.79/13.91 


125/135/84 


0.71 


27.2 


500 


85 


09:45:17 


43.41/14.48/28.93 


82/ 82/ 83 


0.91 


10.3 


500 


108 


10:03:19 


61.40/51.66/13.84 


125/135/85 


0.81 


26.8 


500 


86 


10:19:23 


40.32/13.45/26.87 


86/ 86/ 86 


0.64 


27.5 


500 


109 


10:05:31 


61.18/51.52/13.77 


126/136/85 


0.90 


26.3 


500 


87 


10:21:34 


40.09/13.37/26.72 


86/ 86/ 86 


0.68 


26.1 


500 
















88 


10:23:45 


39.87/13.30/26.57 


86/ 86/ 86 


0.72 


24.6 


500 

















Of Cen A 



2010 Apr 10 


2010 Apr 11 


Cll 


05:57:07 


47.21/15.74/31.46 73/73/73 


0.85 


35.2 


2500 


C16 


05:21:35 


47.66/15.89/31.77 67/67/67 


-99.99 


24.5 


2500 


C12 


06:55:08 


45.81/15.28/30.53 86/86/86 


0.78 


30.6 


2500 


C17 


06:12:28 


46.82/15.61/31.21 77/77/77 


0.67 


20.1 


2500 


C13 


07:46:14 


43.98/14.67/29.32 97/ 97/ 97 


0.67 


21.5 


2500 


C18 


07:00:07 


45.53/15.18/30.35 88/88/88 


0.58 


31.5 


2500 


C14 


09:16:12 


39.78/13.27/26.52 119/118/119 


0.71 


13.8 


2500 














C15 


10:02:42 


37.47/12.49/24.97 131/131/131 


0.69 


13.8 


2500 















Of Cen B 



2010 Apr 10 



2010 Apr 11 



C19 08:58:05 68.07/56.52/13.50 143/150/115 0.54 32.9 2500 
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Fig. C.l. Comparison between the interferometric observables measured in 2009 and those from the reconstructed images within 
the CO line profile shown in Fig.[5l The first, second, third, and fourth columns show the comparison for the continuum, blue wing, 
line center, and red wing, respectively. The filled circles in the top row (a-d) show the positions within the CO line profile. In the 
remaining panels, the observed data and those from the image reconstruction are plotted by the dots and triangles, respectively. The 
reduced ;^^ values including the visibilities, phases, and closure phases, are given in the panels in the bottom row. 
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Fig. C.2. Comparison between the interferometric observables measured in 2010 and those from the reconstructed images within 
the CO line profile shown in Fig. [71 The panels are shown in the same manner as in Fig. lC.ll 



